Seven parton amplitudes from recursion relations 



Daniel de Florian* and Jose Zurita^ 



Departamento de Fisica, FCEYN, Universidad de Buenos Aires, 
(1428) Pabellon 1 Ciudad Universitaria, Capital Federal, Argentina 



Abstract 



We present the first calculation of two-quark and five-gluon tree amplitudes using on-shell 
recursion relations. These amplitudes are needed for tree level 5-jet cross-section and an 
essential ingredient for next-to-leading order 4-jet and next-to-next-to-leading order 3-jet 
production at hadronic colliders. Very compact expressions for all possible helicity configurations 
are provided, allowing for direct implementation in Monte-Carlo codes. 
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1 Introduction 



With the high rate for Standard Model QCD background processes at the forthcoming Large 
Hadronic Collider (LHC), the calculation of multi-jet (particle) production cross-sections becomes 
an essential tool for the discovery of new physics. In order to achieve a good precision within the 
framework of perturbative QCD, next-to-leading order (NLO) calculations are usually needed. 
This is in general a very complicated task in the case of hadronic colliders due to the increasing 
number of involved partons and, for that reason, such a level of accuracy has been reached only 
for a few processes. In several the cases, one still has to rely on tree level calculations. It turns out 
that even at the lowest order the situation starts to be complicated at the level of 5-jets, involving 
the calculation of a few thousand Feynman diagrams. 

After profiting from the great simplifications PUI2] coming out of the combination of the helicity 
method jSl IH and the application of color decomposition rules for amplitudes E] , the task can 
be performed using automatized algorithms, like MadEvent 8J. Even though feasible, that implies 
some non-negligible CPU time for the computation of a few million events, as can be needed for 
simulations. Counting with analytical compact expressions for the amplitudes would certainly be 
a solution for this serious inconvenience. 

Fortunately, the situation has drastically improved during the last couple of years. After the 
pioneering proposal of Witten 9J about the relation between tree level amplitudes and strings 
in twistor space, it was possible to formulate a set of rules to compute gauge amplitudes by 
simple recursion relations involving only "scalar propagators" and the maximally helicity-violating 
(MHV) amplitudes, those where only two particles have a different polarization from the rest. 
Strictly speaking the first set of relations presented by Cachazo, Svrcek and Witten (CSW) [TU] 
involved the off-shell continuation of the MHV amplitudes, situation improved by the proposal 
of Britto, Cachazo and Feng (BCF) ^1] (later confirmed by the same authors and Witten |12j). 
after showing that the usual on-shell MHV amplitudes become the key ingredient when complex 
continuation of some of the external momenta is allowed. The simplicity of the BCFW method 
allows to obtain very compact expressions for those amplitudes, explicitly exposing the high degree 
of symmetry hidden in the framework of direct Feynman diagram calculations. 

The new method, initially considered for pure gluon amplitudes, has been successfully extended 
to account for the presence of massless quarks O E] , Higgs boson [T3] , massive gauge bosons 
photons ^7] and even massive fermions |18| . In the case of pure gluonic processes, helicity 
amplitudes involving up to eight particles have been computed. For those involving also massless 
fermions calculations have been performed up to six particles |19| I2U|. Furthermore, recent progress 
has been done to extend the validity of the recursion relations to one-loop amplitudes in QCD 

HDI221l2ni. 

In this paper, we present the first calculation of the full set of helicity amplitudes involving 
a quark- ant iquark pair and five-gluons, needed for the computation of tree-level five-jet cross- 
sections, and an ingredient for the real part of NLO(NNLO) results for four(three)-jet observables 
in hadronic collisions. 

This paper is organized as follows: in Section 2 we review the main ingredients of the BCFW 
formulation and recall its limitations when fermions are present. In Section 3 we introduce the 
main results for the complete set of qq5g helicity amplitudes, while in Section 4 we present our 



1 



conclusions. 



2 Color decomposition, helicity and BCFW 

The color decomposition for a qq pair and n gluons at tree level is particularly simple. The 
amplitude Mn can be written in terms of the partial amplitude An as PQ 

Mi°\k l ,\,a i )=g n - 2 (^ aCT(3) ---^0|4° ) (lJ 1 ! 25 2 !( T(3 A3 ) ! ...,a(n A -)), (1) 

where S" n _2 is the group of permutations of n — 2 symbols, with 1 representing the quark with 
color i\ and 2 the antiquark with color j'2. The upper- index A; indicates the helicity of particle I 
carrying momentum fcj. The normalization for the color matrices in the fundamental representation 
is Tr (T a T b ) = S ab . 

In the framework of the helicity formalism jHlElE], with the spinors denoted as 

= \kf) = Mh) (^| = (kt\ = Mklj, (2) 

the partial amplitudes can be written in terms of the spinors inner-products 

and a few simple combinations of them, like 

(i\Pa\j] = (ia)[aj] 
(i\PaPb\j) = (ia)[ab)(bj) . (4) 

In our convention all particles are considered to be outgoing and, following the QCD literature 



[UEj, we fix the sign of the inner products such that = Sj- 



The BCFW recurrence relation is based on the analytical properties of the amplitude when 
the spinors of two external legs (denoted by j and /) are shifted as 



Ij> = li> 
\j] = \j]-z\i) 

\t) = \l)+z\j) 

\i] = \i] . 

After this shift, the BCFW formula simply reads 

r,s A=± 

' A^ s+r (K~ x , (s + 1) W), . . . , l\ . . . , (r - , 



(5) 



(6) 



K 2 



* When comparing with results obtained using the string-like conventions just notice that [ij] carries the opposite 
sign 
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where K rs = k r + . . . + kj . . . + k s and the (complex) shift variable z takes the value 

(j\K rs \l}- 



(7) 



At this point we should make a few remarks on Eq.([fi|). First of all, each term is the product 
of two helicity amplitudes with a fewer number of particles and a propagator. The sum over r 
and s is not actually a sum over all of their possible values, but instead over all the possible 
configurations where the j particle belongs to one of the amplitudes and the I particle to the other 
one. For future reference, we shall call the amplitude with the j particle "upper amplitude" and 
the one including the I particle "lower amplitude" . There are certain restrictions to the massless 
particles that can be used as the reference lines j and I; in general they can not be chosen as 
(Aj,A;) = (+,—). Furthermore, quarks and antiquarks of the same flavor can not be chosen if 
they are adjacent and for adjacent quarks and gluons the helicities should better be opposite [T9"] . 

We should also notice that the sum includes amplitudes involving only 3 on-shell partons. 
Because of helicity conservation these amplitudes would vanish if the momenta were not shifted. 
It is straightforward to show that, after the shift in Eq.(j5J), only the 3 parton upper MHV and the 
3-parton lower MHV amplitudes become non-zero. Therefore, the ggg and qqg MHV amplitudes 
which with our phase conventions read 

Az (12)(31) 

4 0) (l, + ,2 5 ,3 9 ) = ^, (8) 



are the key ingredients of the recursion relations. The corresponding MHV amplitudes can be 
obtained from those above by flipping the helicities applying parity inversion ((ij) — > [ji] and an 
extra factor of -1 for each pair of quarks participating) and charge conjugation plus reflection and 
cyclic symmetries of the amplitudes. Using the recursion relations it is possible to construct the 
tree level amplitude for n-partons just by conveniently iterating these building blocks. 



3 qq5g helicity amplitudes 

There are in principle 2 7 different helicity amplitudes for this process but half of them, those with 
the quark and the antiquark carrying the same helicity, are trivially vanishing for massless particles. 
Furthermore, it is enough to present the results for one of the two possible combinations of qq 
helicities (we choose here lj" 2^); the other can be obtained by parity and charge conjugation. Out 
of the remaining 32 amplitudes, those with all gluons with the same helicity (two) are vanishing 
and other 10 are either MHV or MHV amplitudes [23], simply reading 

(o) f + + +) _ (2t) 3 (U) 

■"-1 \ L q i i °g > • • • 5 b g • • ■ ; ' g J ~ y-r7 



^7 2g j 3 g , . . . , . . ■ , 7 g ) — i _ [7 rj ; ■ — , (9) 
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where % represents the gluon with the opposite helicity to the others. 

Therefore, there are only 10 non-trivial NMHV, corresponding to three gluons with helicity plus 
and two with helicity minus, amplitudes to be computed. Again, the 10 NMHV amplitudes can be 
obtained by the discrete symmetries P and C. Further simplifications in the number of independent 
amplitudes could be achieved by applying supersymmetric relations. We rather present the explicit 
results for those 10 amplitudes in order to provide the most compact expressions for direct use. 

The use of the BCFW formula for qq5g amplitudes involves the appearance of, at most, four 
different arrangements in the recursion. This is shown in Fig.(£Q) for the (1+, 2q , 3 + , 4 + , 5 + , 6~, 7~) 
helicity configuration, where gluons 6 and 7 were chosen as lines j and I in Eq.(J3J, respectively. 
As known, a smart election for the reference lines can result into more compact expressions for 
the amplitudes. 




Figure 1: Diagrams contributing to A\' (1+, 2 5 , 3 + , 4 + , 5 + , 6 ,7 ) 



In this case diagram (a), denoted as (2, 3, 4, 5, 6|7, 1) vanishes because the two fermions appear 
in different subamplitudes, fixing the helicity of the propagator and selecting a 3-parton lower MHV 
amplitude. The remaining diagrams ( (3, 4, 5, 6|7, 1, 2), (4, 5, 6|7, 1, 2, 3), and (5, 6|7, 1, 2, 3, 4)) are 
simply products of MHV and/or MHV amplitudes, so each of them contributes with a single term. 
Our final result reads 



,(o) fl+ 2 - 3+ 4+ ,+ 6 _ = (6|7 + 2|1] 

7l «'<?' ' ' ' ' } ~ s n2 [71][12](34)(45)(56)(3|l + 2|7] 



(2|6 + 7|5] 3 (2|(l + 3)(4 + 5)|6) 



(10) 



s 567 {12> (23) (34) [56] [67] (4|5 + 6|7] s 123 s 4 56 (12) (23) (45) (56) (3| 1 + 2|7](4|5 + 6|7] ' 

where Sijk = (p% + Pj + Pk) 2 - Factorization properties of the amplitudes in the collinear and soft 
limits provide stringent consistency checks to them. For example, when gluon 6 is soft, Eq. (jlO|) 
becomes 

4 0) (1+ 2" 3 + ,4 + ,5 + ,6-,7-) ^ ( , w "f 27 ^ w ^ (- J^L), (11) 
7Kql q1 ' ' ' ' ; V( 12 ) (23) (34) (45) (57)7 V [56][61] 7 ' V 1 

i.e., just the product of the six particle amplitude Af\l+ , 2^, 3 + , 4 + , 5 + , 7 _ ) times the eikonal 
factor for the emission of a soft gluon with negative helicity. 

In our search for compact expressions, we can show that the NMHV 7-parton amplitudes have 
at most six terms. In any amplitude, diagrams like (a) contribute with (at most) three terms 
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(since it involves a 6 parton NMHV amplitude QjJ and each of the other diagrams adds a 
single term, if they don't vanish. 



The amplitude corresponding to the ordering + H 1 — on the gluon helicities is computed 

choosing j = 3 and I = 4 as reference lines, resulting 

4 0) (1+ 2r,3 + ,4+,5-,6 + ,7-) = 

(5|3 + 4|1] 3 (57) 4 

[12](34)(45)(56)(67)(5|6 + 7|1](3|(1 + 2) (6 + 7)|5)(7|(1 + 2) (3 + 4)|5) 
(27) 3 (5|3 + 4|6] 4 

+ S712S3 4 5 (12)(34)(45)(3|4 + 5|6](2|7+l|6](7|(l + 2)(3 + 4)|5) 

(25) 3 [16] 3 (5|7+1|6] (12) 

s 671 (23) (34) (45) [67] [71] (2 1 7 + 1 1 6] (5 1 6 + 7 1 1] 
(27) 3 [46] 4 



+ 



S456 (12)(23)[45][56](3|4 + 5|6](7|6 + 5|4] 

(2|1 + 3|4] 3 (57) 4 

si23S 567 (12)(23)(56)(67)(7|6 + 5|4](3|(1 + 2)(6 + 7)|5) ' 



where we have splitted the result in the following order: the first three terms come from (6, 7, 1, 2, 3|4, 5), 
the next single term from (7, 1, 2, 3|4, 5, 6) and the last term from (1, 2, 3|4, 5, 6, 7). The contri- 
bution from (2, 3|4, 5, 6, 7, 1) vanishes because the "fermion propagator" selects an upper MHV 
amplitude. 

The amplitude corresponding to the ordering -H h on the gluon helicities is also computed 

choosing j = 3 and / = 4 as reference lines, obtaining 



+ 



(o) /a 

[17] 2 (25) 3 (5|1 + 6|7] 



^(1+2-3+4+5-6-7 



+ 
+ 



s 671 (23) (34) (45) [67] (5 1 6 + 7 1 1] (2 1 7 + 1 1 6] 

(2|(7+l)(3 + 4)|5) 3 (l|(2 + 7)(3 + 4)|5) 

s 712S3 45(71)(12)(34)(45)(3|4 + 5|6](2|7 + l|6](7|(l + 2)(3 + 4)|5) 

(5|4 + 3|1] 3 (56) 3 (13) 

[12](34)(45)(67)(5|6 + 7|1](3|(1 + 2)(6 + 7)|5)(7|(1 + 2)(3 + 4)|5> 

(2|5 + 6|4] 3 (1|5 + 6|4] 
S 45 6 (71)(12)(23)[45][56](3|4+5|6](7|5 + 6|4] 
(2|1 + 3|4] 3 (56) 3 



s 567 s 123 (12)(23)(67)(7|5 + 6|4](3|(1 + 2)(6 + 7)|5) 



The first three terms come from (6, 7, 1, 2, 3|4, 5), the fourth from (7, 1, 2, 3|4, 5, 6) and the fifth from 
(1, 2, 3|4, 5, 6, 7). Again, because of the "fermion propagator" and the helicities of the particles 
involved, (2, 3|4, 5, 6, 7, 1) has a null contribution to this result. 
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The amplitude corresponding to the ordering H h H — is obtained selecting j — 5 and / = 6 

as reference lines, 

4 0) (l+,2 9 r,3+,4-,5 + ,6 + ,7-) = 
(24) 3 (7|5 + 6|l] 3 (4|(2 + 3)(5 + 6)|7) 

s 2 34W23)(34>(56>(67><4|2 + 3|1](5|6 + 7|1](2|(3 + 4)(5 + 6)|7) 

(27) 3 (7|5 + 6|3] 4 

" (12)[34](56)(67)(7|1 + 2|3](7|5 + 6|4](5|(3 + 4)(1 + 2)|7)(2|(3 + 4) (5 + 6)|7) 

[13] 3 (47) 4 

~ s 123 [12](45)(56>(67>(4|2 + 3|1](7|1 + 2|3] (14) 

(24) 3 [16] 3 (4|7+ 1|6] 
~ s 671 (23)(34)(45)[67][71](2|7+l|6](5|6 + 7|l] 

(27) 3 (4|5 + 3|6] 4 

~ s 712 S34 5 <12)(34>(45)(3|4 + 5|6](2|7+l|6](5|(3 + 4)(l + 2)|7) 
(27) 3 [56] 3 

s 4 5 6 (12)(23)[45](3|4 + 5|6](7|6 + 5|4] " 

This formula was written in the following order: the first three terms come from (1, 2, 3, 4, 5|6, 7), 
the next one from (2, 3, 4, 5|6, 7, 1), a single term from (3, 4, 5|6, 7, 1, 2) and the last one comes 
from (4,5|6,7, 1,2,3). 

The ordering H 1 h can be obtaining setting j = 6 and I = 7, 

4 0) (l+,2r,3+,4-,5 + ,6-,7 + ) = 

(16)(24) 3 (4|2 + 3|5](6|7+1|5] 3 

s 6 nS234(23)(34)(67)(71)(l|6 + 7|5](2|3 + 4|5](6|(7+l)(2 + 3)|4) 

(16)(26) 3 [35] 4 

+ s 3 45(67)(71)(12)[34][45](6|4 + 5|3](2|3 + 4|5] 

(16)(46) 4 (6|7+1|3] 3 

+ (45)(56)(67)(71)(6|4 + 5|3](6|7 + 1|2](6|(4 + 5) (2 + 3)|1)(6|(7 + 1)(2 + 3)|4) (15) 

[27][17] 2 (46) 4 

+ s 712 [12](34)(45)(56)(3|l + 2|7](6|7 + 1|2] 

(1|2 + 3|7](2|1 + 3|7] 3 (46) 4 

+ Sl23 S456(12)(23)(45)(56)(4|5 + 6|7](3|1 + 2|7](6|(4 + 5) (2 + 3)|1) 
, (14)(24) 3 [57] 4 

s 5 67(12)(23)(34)[56][67](4|5 + 6|7](1|6 + 7|5] ' 

where the first three terms come from (2, 3, 4, 5, 6|7, 1), the fourth term from (3, 4, 5, 6|7, 1, 2), the 
fifth from (4, 5, 6|7, 1, 2, 3) and the sixth term from (5, 6|7, 1, 2, 3, 4). 
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The amplitude corresponding to the ordering H 1 — h is obtained by choosing j = 6 and 

I = 7 as reference lines, 



4 0) (l+,2r,3 + ,4-,5-,6+,7+) = 

, (4|(2 + 3)(6 + 7)|l)( S671 ) 2 (24)3 

+ S234 (23)(34)(67)(71)(l|6 + 7|5](2|3 + 4|5](4|(2 + 3) (7 + 1)|6) 

(1|4 + 5|3](2|4 + 5|3] 3 

+ S34 5 (67) (71) (12) [34] [45] (2|3 + 4|5](6|5 + 4|3] 

(16)(45) 3 (6|7+1|3] 3 

~ (56)(67)(71)(6|7+ 1|2](6|5 + 4|3](6|(4 + 5)(2 + 3)|1)(4|(2 + 3)(7 + 1)|6) (16) 

, [27][17] 2 (45) 3 

S712 [12] <34> <56> <3 ] 1 + 2 1 7] <6 1 7 + 1 1 2] 

(1|2 + 3|7](2|1 + 3|7] 3 (45) 3 

+ Sl23S4 56(12)(23)(56)(3|l + 2|7](4|5 + 6|7](6|(4 + 5) (2 + 3)|1) 

(14)(24) 3 [67] 3 

+ S5 67 (12)(23)(34)[56](4|5 + 6|7](l|6 + 7|5] ' 

where we have written the result in the following order: the first three terms come from (2, 3, 4, 5, 6|7, 1), 
a single term from (3, 4, 5, 6|7, 1, 2), the next term from (4,5,6)7,1,2,3) and the last term from 
(5,6|7,1,2,3,4). 

The ordering — h + H — , using j = 7 and 1 = 1, results 
4 0) (l+,2r,3-,4+,5 + ,6 + ,7-) = 

(7|5 + 6|4] 3 (27) 3 

(12)[34](56)(67)(7|1 + 2|3](2|(3 + 4)(5 + 6)|7)(7|(1 + 2) (3 + 4)|5) 

(3|4 + 5|6] 3 (27) 3 

+ s 712 S345 < 12) <34> <45> <2| 7 + 1 1 6] <7| (1 + 2) (3 + 4)|5) 

(7|1 + 3|2](7|2 + 3|1] 3 (17) 

s 123 [12] [23] (45) (56) (67) (7| 1 + 2|3] (4|3 + 2| 1] 
(7|(5 + 6)(2 + 4)|3)(7|5 + 6|l] 3 (23) 3 

s 2 34S 5 67(23)(34)(56)(67)(4|3 + 2|l](5|6 + 7|l](7|(5 + 6)(3 + 4)|2) 
(3|7+1|6](23) 2 [16] 3 

s 671 (34) (45) [67] [71] (2|7 + 1 16] (5|6 + 7| 1] ' 

The different contributions appear in the following order: the first two terms come from (3, 4, 5, 6, 7|1, 2), 
the third term comes from (4, 5, 6, 7|1, 2, 3), the fourth term from (5, 6, 7|1, 2, 3, 4) and the last term 
from (6,7|i,2,3,4,5). 
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The amplitude corresponding to the ordering — I 1 — h is obtained by selecting j = 5 and 

I = 6 as reference lines, 

4 0) (l+ ! 2r,3-,4 + ,5-,6 + ,7 + ) = 

(15)[24](5|2 + 3|4] 3 



+ 
+ 
+ 
+ 



s 23 4[23][34](56)(67)(71)(5|3 + 4|2](1|2 + 3|4] 

(5|6 + 7|2](5|6 + 7|1] 2 (35) 4 

[12](34)(45)(56)(67)(5|3 + 4|2](3|(1 + 2)(6 + 7)|5)(5|(3 + 4)(1 + 2)|7) 

(13)(23) 2 (5|6 + 7|4] 4 (18) 

s 5 67Si2 3 (12)(56)(67)(7|6 + 5|4](1|2 + 3|4](3|(1 + 2)(6 + 7)|5) 

(1|7 + 2|6](2|7+ 1|6] 2 (35) 4 

s 712 s 3 4 5 (71)(12)(34)(45)(3|4 + 5|6](5|(3 + 4)(1 + 2)|7) 

(13)(23) 2 [46] 4 

5456(71) (12) [45] [56] (3|4 + 5|6](7|5 + 6|4] ' 



where we have written the result in the following order: the three terms from (1, 2, 3, 4, 5|6, 7), a 
single term from (3, 4, 5|6, 7, 1, 2) and the single term from (4, 5|6, 7, 1, 2, 3). The contribution from 
(2, 3, 4, 5|6, 7, 1) vanishes because of the "fermion propagator" and the helicities of the particles 
involved. 

The ordering — h H h is solved choosing j = 3 and I = 4 

4 0) (l+,2r,3-,4 + ,5 + ,6-,7 + ) = 

(16)(3|4 + 5|2](6|(7+1)(4 + 5)|3) 3 



+ 
+ 



S345S 67 i(34)(45)(67)(71)(5|4 + 3|2](6|7+l|2](l|(6 + 7)(4 + 5)|3) 

[27][17] 2 (36) 4 

+ s 712 [12] (34) (45) (56) (6 1 7 + 1 1 2] (3 1 1 + 2 1 7] 

(13)(23) 2 (3|4 + 5|7] 4 

~~ (12)(34)(45)[67](3|4 + 5|6](3|1 + 2|7](3|(1 + 2)(6 + 7)|5)(1|(6 + 7) (4 + 5)|3) (19) 

(13)(23) 2 [45] 3 

S 4 56 ( 71)(12)[56](3|4 + 5|6](7|5 + 6|4] 

(13)(23) 2 (6|5 + 7|4] 4 

S567Si 23 (12)(56)(67)(l|2 + 3|4](7]5 + 6|4](3|(1 + 2)(6 + 7)|5) 
(16)[24](6|2 + 3|4] 3 
~ s 234 [23][34](56)(67)(71)(l|2 + 3|4](5|4 + 3|2] " 

This expression was splitted in the following order: the first three terms come from (6, 7, 1, 2, 3 1 4, 5), 
the fourth from (7, 1, 2, 3|4, 5, 6), the fifth term from (1, 2, 3|4, 5, 6, 7) and the sixth term from 
(2,3|4,5,6,7,1). 

Finally, the non-alternating amplitude h + + (from j = 4 and I = 5) turns out to be 

4 0) (1 + ,2" 3-,4-,5 + ,6 + ,7 + ) = 



(l|(2 + 7)(5 + 6)|4)(2|(7+l)(5 + 6)|4) 2 
' ' ' S456S 7 i2(71)(12)(45)(56)(6|5 + 4|3](7|l + 2|3] 



(4|1 + 3|2](4|2 + 3|1] 2 (1|3 + 4|5](2|3 + 4|5] 2 



(20) 



s 123 [12][23](45)(56)(67)(7|l + 2|3] s 345 <67><71)<12>[34][45]<6|5 + 4|3] 
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where we have written the result in the following order: the first two terms come from (7, 1, 2, 3, 4|5, 6), 
and the last one from (3, 4|5, 6, 7, 1, 2). Contributions from (1, 2, 3, 4|5, 6, 7) and (2, 3, 4|5, 6, 7, 1) 
vanish because of helicity conservation. 

As a byproduct of this calculation, several QED amplitudes can be obtained form our re- 
sults, by turning any number of gluons into photons. In order to disguise gluons as photons one 
has to replace in Eq.(JTJ) the SU(3) color group of QCD by the U(l) group of QED. The color 
decomposition for two quarks, m photons and r gluons [23] becomes 

M^{h, Aj, at) = g r (V2eQ q ) m ]T (T a ^ . . . TM^))| 

aeSr (21) 

^(1^,2^,^(3^), . . .,a g ((r + 2) x <^), 7r+3 , . . . , 7n ) , 

where S r is the group of permutations of r symbols with 1 being the quark with color i±, 2 the 
antiquark with color j'2, and A/ represents the helicity of the / particle carrying momentum ki. This 
formula can be easily obtained from Eq. by converting the last m gluons into photons. In order 
to make this conversion, one has to replace the factor gT a with ^/2eQ g I, where / is the identity 
matrix. Since / commutes with the T a matrices, all possible photon-quark couplings contribute to 
the same color structure. Therefore to obtain the amplitudes involving photons, one has to sum 
over all the permutations where the photon "moves" in the amplitude while the gluons remain 
fixed. 

From the previous results one can change any number of gluons into photons, using the proce- 
dure described in the previous paragraph, allowing to obtain the following QED/QCD amplitudes: 
qq4g>y, qq3g2^, qq2g3>y, qqg4>y and qq5>y . 



4 Conclusions 

In this paper we presented all seven parton tree level NMHV amplitudes involving a fermionic 
pair and five gluons, obtained by use of the BCFW recursion relations. With the knowledge of 
these amplitudes the full set of helicity amplitudes for the two quarks plus five gluons process is 
available. The trivial MHV amplitudes are giving by the Parke- Taylor formulae and the NNMHV 
amplitudes are also NMHV so they can be obtained by performing the adequate combination of 
P and C discrete symmetries over our results. We should emphasize that the results presented in 
this paper have been checked in all possible collinear and soft limits, setting and stringent test for 
the correctness of the amplitudes 

Furthermore, by making simple replacements in the color decomposition formula, one can 
obtain several seven parton QED/QCD amplitudes involving two quarks, m photons and r gluons 
(where m + r = 5). 

These amplitudes are a main ingredient for the calculation of multijets cross sections in 
hadronic colliders. As expected, we have obtained very compact expressions for the amplitudes, 
allowing for a more convenient implementation in computer codes than those obtained from au- 
tomatic tree level computations methods. 

■'"Stricktly speaking they are correct up to terms that must vanish in all possible soft and collinear limits, which 
are very unlikely to exist. 
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